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Understanding  the  physics  and  predicting  transition  from  laminar  to 
turbulent  flow  remains  one  of  the  central  problems  in  fluid  dynamics. 
Theoretical  and  experimental  studies  have  shown  that,  if  free  stream 
disturbances,  wall  vibrations,  etc.,  are  sufficiently  small,  the  appearance 
and  growth  of  Tollmien-Schlichting  waves  can  dominate  the  process  of  boundary 
layer  transition.  Recent  theoretical  studies  suggest  that  the  Tollmien- 
Schlichting  waves  do  not  constitute  a complete  set  of  solutions  of  the 
Orr- Sommer f eld  equation  and  hence  cannot  be  used  to  represent  an  arbitrary 
disturbance.  In  addition,  it  has  long  been  recognized  that  the  Tollmien- 
Schlichting  waves  are  largely  confined  to  a region  within  two  or  three 
boundary  layer  thicknesses  of  the  solid  wall.  The  excitation  of  these 
waves  by  free  stream  disturbances  which  do  not  penetrate  the  boundary 
layer  is  an  unsolved  problem. 

In  this  paper  it  is  shown  that  the  Orr-Sommerfeld  equation  has  a 
continuous  spectrum,  in  addition  to  the  discrete  spectrum  of  the  Tollmien- 
Schlichting  waves,  for  any  flow  in  an  unbounded  domain  such  as  a boundary 
layer,  jet  or  wake.  Explicit  formulae  for  the  location  of  the  continuum 
in  the  complex  wave  speed  plane  are  given.  The  continuous  spectrum  and 
the  continuum  eigenfunctions  are  calculated  for  two  sample  problems:  the 

Blasius  boundary  layer  and  the  two-dimensional  laminar  jet.  It  is  shown 
that  these  modes  arc  free  stream  modes;  far  from  the  boiuidary  or  jet  or 
wake  axis  they  are  vorticity  waves.  Three  mechanisms  are  suggested  by 
which  free  stream  vorticity  could  cause  transition  in  a shear  flow  by 
exciting  the  continuum  modes.  These  mechanisms  are  (1)  quasi-steady 
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distortion  of  the  mc-in  profile;  (2)  time- dependent  instability  of  the  mean  j 
flow  plus  a continuum  mode;  and  (3)  the  effects  of  nonparallel  flow  correc-  j 
tions  to  the  continuum  spectrum  and  eigenfunctions.  t 


Further  research  is  needed  to  show  how  a "patch"  of  vorticity  in  the 
free  stream  (a  wave  pocket)  can  be  represented  in  terms  of  the  continuvun 
eigenfunctions;  to  elucidate  the  connection  between  these  modes  and  the 
forced  disturbances  studied  by  others;  and  to  examine  the  postulated 
transition  mechanisms  listed  above.  We  have  begun  a theoretical  study 
of  these  problems.  It  would  be  desirable  if  experimentalists  would  also 
look  for  these  modes. 


1.  INTRODUCTION 


The  study  of  the  stability  of  laminar  flow  against  externally  imposed 
disturbances  and  the  connection  of  such  instability,  if  it  exists,  with  the 
transition  from  laminar  to  turbulent  flow  has  been  studied  for  nearly  a 
century  and  remains  one  of  the  central  problems  in  fluid  ctynamics.  Ever 
since  Reynolds'  classic  experiments  (1883)  it  has  been  conjectured  that  tran> 
sition  from  laminar  to  turbulent  flow  is  the  result  of  an  instability  in  the 
laminar  flow.  Beginning  with  Rayleigh  (1880),  there  have  been  many 
theoretical  studies  which  have  attempted  to  predict  under  what  conditions 
small  disturbances  in  the  velocity  profile  would  grow  or  decay. 

The  main  interest  in  stability  calculations  centers  on  the  stability 
of  boundary  layer  flows  and  related  flows  such  as  jets,  wakes,  free  shear 
flows;  Couette  and  particularly  Poiseuille  flow  have  had  the  status  of 
"canonical"  problems.  Tollmien  (1931),  Schlichting  (1951)  and  Lin  (1955) 
made  notable  advances  in  the  analytical  calculation  of  the  stability  of  these 
parallel  shear  flows  to  two-dimensional  wavelike  disturbances,  the  Tollmien- 
Schlichting  waves  (hereafter  called  T-S  waves).  The  restriction  to  two-dimensional 
waves  was  justified  by  a result  of  Squire  (1933) , who  showed  that  if  an  oblique 
wave  is  unstable  at  some  Reynolds  number,  then  a two-dimensional  T-S  wave  is 
unstable  at  a smaller  Reynolds  number.  Thus,  in  order  to  find  the  critical 
Reynolds  number  it  suffices  to  consider  only  two-dimensional  T-S  waves. 

Although  the  theory  was  highly  developed,  it  was,  prior  to  the  early  1940* s, 
generally  regarded  as  irrelevant  to  the  actual  mechanics  of  transition  because 
T-S  waves  had  never  been  observed  experimentally,  and  it  appeared  that  any 
disturbance,  whether  in  the  free  stream,  in  the  boundary  layer,  or  on  the 
boundary,  could  catise  a transition  if  it  were  large  enough  or  if  the  Reynolds 
number  were  high  enough. 


The  experimental  studies  of  Schubauer  and  Skramstad  (1941,  1947,  1948) 
first  showed  that,  if  the  free  stream  disturbances,  wall  vibrations,  etc., 
were  sufficiently  small,  the  appearance  and  growth  of  T-S  waves  could  play  a 


dominant  role  in  boundary  layer  transition.  Schubauer  and  Skramstad  found 
generally  good  agreement  between  the  predictions  of  linear  stability  theory 
and  the  experimental  results,  not  only  for  the  shape  of  the  neutral  stability 
curve  and  the  phase  speeds  of  the  T-S  waves,  but  also  for  the  shape  of  the  T-S 
waves  as  a function  of  distance  from  the  boundary. 

While  the  work  of  Schubauer  and  Skramstad  showed  the  essential  correct- 
ness of  linear  stability  theory,  the  calculat_ons  of  the  detailed  stability 
characteristics  for  a particular  profile  remained  quite  delicate  in  that  the 
results  obtained  were  sensitive  to  the  detailed  shape  of  the  profiles  and  the 
approximations  used  in  the  asymptotic  analysis.  It  is  only  the  development  of 
sophisticated  numerical  methods  and  the  high  speed  computers  to  implement  them 
(Kaplan,  1964;  Gallagher  and  Mercer,  1964;  Landahl,  1966;  Radbill  and  Van 
Driest,  1966;  Wazzan,  Okamura  and  Smith,  1968;  Lessen  et  al.,  1968;  Grosch  and 
Salwen,  1969;  Orszag,  1971;  Salwen  and  Grosch,  1972)  that  has  made  these 
calculations  relatively  routine. 

It  is  usually  assumed  that  the  eigenvalue  problem  of  linear  stability 
theory  has  an  infinite  set  of  discrete,  eigenvalues  and  a corresponding 
infinite  set  of  eigenfunctions.  Assuming  the  completeness  of  this  set  of 
eigenfunctions,  an  arbitrary,  initial  disturbance  which  satisfies  the 
boundary  conditions  can  be  expanded  in  terms  of  them.  We  know  of  no  general 
proof  of  this  completeness  ''theorem." 

Almost  ail  stability  calculations  have  been  concerned  with  finding  the 
first  or  least  stable  eigenvalue.  It  should  be  noted  that,  if  expansion 
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(spectral)  techniques  (Grosch  and  Salwen,  1969;  Orszag,  1971;  Salwen  and 
GroSch,  1972)  are  used  to  solve  the  eigenvalue  problem  in  a finite  domain, 
a number  of  the  higher  eigenvalues  are  found  numerically  along  with  the 
first.  If  a sufficient  number  of  terms  is  used  in  the  expansion,  these 
numerically  calculated  eigenvalues  will  overlap  the  asymptotic  (high  order) 
distribution  of  eigenvalues,  for  which  analytic  formulae  may  be  derived, 
and  the  complete  spectrum  may  be  found.  Apart  from  any  difficulties  in  the 
analysis  or  calculation,  this  concentration  on  the  first  eigenvalue  may  be 
due  to  the  facts  that  1)  experiments  suggest  that  there  is  only  a single 
unstable  mode  and  2)  for  those  flows  for  which  the  higher  modes  have  been 
calculated,  plane  Couette  (Gallagher  and  Mercer,  1964),  plane  Poiseuille 
(Grosch  and  Salwen,  1969;  Orszag,  1971),  pipe  Poiseuille  (Lessen  et  al., 

1968;  Salwen  and  Grosch,  1972) , these  higher  modes  have  been  found  to  be 
highly  damped. 

Very  recently  Jordinson  (1971)  and  Mack  (1975),  using  different 
numerical  methods,  have  calculated  the  higher  eigenvalues  of  the  Orr- 
Sommerfeld  equation  for  Blasius  flow.  Jordinson  calculated  eigenvalues  for 
the  spatially  and  temperally  growing  or  decaying  modes  for  a single  Reynolds 
number,  and  a single  wave  number  (temporal  problem)  or  a single  frequency 
(spatial  problem).  Mack  calculated  the  eigenvalues  for  a number  of 
different  values  of  the  wave  number  and  Reynolds  number  for  the  tenqjoral 
stability  problem.  While  there  is  some  disagreement  in  the  number  and 
location  of  the  eigenvalues  found  by  Jordinson  and  Mack,  they  both  agree  that 
at  any  Reynolds  number  there  are  only  a finite,  and  small,  number  of  discrete 
eigenvalues.  A finite  set  of  eigenfunctions  cannot  be  a complete  set.  How, 
then,  can  an  arbitrary  initial  disturbance  be  expanded?  Both  Jordinson  and 
Mack  suggest  that  there  is  a continuous  spectrum. 


We  know  of  no  general  results  concerning  the  spectrum  of  the  Orr- 
Sommerfeld  equation  for  an  arbitrary  parallel  shear  flow  in  an  unbounded 
domain.  For  the  Orr-Sommerfeld  equation  in  a bounded  domain  or  in  the  case 
of  the  inviscid  stability  problem,  for  which  the  governing  equation  is  the 
Rayleigh  equation,  a number  of  specific  problems  have  been  studied  in 
detail,  and  some  general  results  concerning  the  spectrum  are  available. 

Much  is  known  about  the  spectrum  of  the  inviscid  stability  problem 
(the  Rayleigh  problem).  Howard's  circle  theorem  (1961)  states  that  all 
eigenvalues  of  the  Rayleigh  problem  lie  in  a circle  in  the  complex  plane 
whose  diameter  is  the  range  of  undisturbed  velocity.  Case  (1960)  has  studied 
the  spectrum  of  inviscid  plane  Couette  flow  and  shown  that  there  are  no 
discrete  eigenvalues;  there  is  only  a continuum.  In  the  seme  paper  Case 
generalized  this  result  and  showed  that  for  any  mean  velocity  profile  in  a 
finite  domain,  the  Rayleigh  problem  always  has  a continuous  spectrum  in 
addition  to  the  discrete  spectrum,  if  one  exists.  Rosencrans  and  Sattinger 
(1966)  and  Sattinger  (1967)  have  proved,  under  fairly  general  conditions  on 
the  velocity  profile  of  the  mean  flow,  that  the  Rayleigh  problem  in  a finite 
domain  has  only  a finite  number  of  discrete  eigenvalues.  Plane  Couette  and 
Poiseuille  flows  and,  it  appears,  most  other  physically  reasonable  flows  in 
a finite  domain  satisfy  the  conditions  of  this  theorem. 

Case  (1961)  and  Lin  (1961)  have  examined  the  connection  between  the 
spectrum  of  the  invis'^ld  problem  and  that  of  the  viscous  problem  for  a finite 
domain.  Case  has  shoim  that,  if  the  initial  conditions  are  independent  of 
the  viscosity  (Reynolds  number) , the  solution  of  the  viscous  problem 
approaches  the  solution  of  the  inviscid  problem  as  the  viscosity  goes  to  zero 
and,  hence,  has  a continuous  spectrum.  Lin  has  pointed  out  that,  for  the 
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viscous  problem,  the  initial  conditions  depend  in  general  on  the  viscosity  and, 
therefore.  Case's  theorem  does  not  apply.  Lin  has  further  shown  that  the 
Orr-Sommerfeld  equation  has  dnly  discrete  eigenvalues  in  a bounded  domain. 

It  has  been  proved  that  the  Orr-Sommerfeld  equation  has  a complete 
set  of  eigenfunctions  for  the  specific  cases  of  plane  Couette  flow  (Haupt, 

1912)  and  plane  Poiseuille  flow  (Schensted,  1960).  Di  Prima  and  Habelter 
(1969)  have  proved  a completeness  theorem  for  a class  of  nonself-adjoint 
eigenvalue  problems  in  a bounded  domain.  Using  this  theorem  they  showed  that 
the  Benard  problem  (linear  stability  of  a layer  of  fluid  heated  from  below), 
the  Taylor  problem  (linear  stability  of  the  flow  between  rotating  concentric 
cylinders) , and  the  Orr-Sommerfeld  equation  for  any  flow  in  a bounded 
domain  have  a complete  set  of  eigenfunctions. 

None  of  these  results  apply  to  the  Orr-Sommerfeld  equation  if  the  domain 
is  infinite  as  it  is  for,  say,  Blasius  flow.  The  possibilities  for  the 
spectrum  of  the  Orr-Sommerfeld  equation  in  an  infinite  domain  are:  (a)  there 
is  an  infinite  set  of  discrete  eigenvalues  without  a continuum;  (b)  there 
is  an  infinite  set  of  discrete  eigenvalues  with  a continuum;  (c)  there  is 
a finite  number  of  eigenvalues  without  a continuum;  or  (d)  there  is  a finite 
number  of  discrete  eigenvalues  with  a continuum. 

In  this  paper  we  will  not  be  concerned  with  whether  or  not  the  discrete 
spectrum  of  the  Orr-Sommerfeld  equation  is  finite  in  an  infinite  domain.  We 
will  consider  the  existence  of  the  continuous  spectrum  of  the  Orr-Sommerfeld 
equation  for  both  the  temporal  and  spatial  problem.  In  section  2,  it  is  shown 
that  the  Orr-Sommerfeld  equation  has  a continuous  spectrum  for  a wide 
class  of  shear  flows  in  an  infinite  domain.  Formulae  for  the  continuous 
spectrum  in  terras  of  the  wave  number  or  frequency  and  the  Reynolds  number 


are  given  for  these  flows.  Detailed  calculations  of  the  continuous  spectrum 
and  examples  of  the  continuum  eigenfunctions  for  two  specific  flows,  the  Blasius 
boundary  layer  and  the  two-dimensional  laminar  jet,  are  given  in  section  3. 

Also  included  in  section  3 is  a physical  interpretation  of  the  modes.  Section 
4 contains  some  speculation  on  mechanisms  by  which  these  modes  could 
lead  to  transition.  Finally,  section  5 is  a brief  summary  of  our  results. 

In  a later  paper  (Part  II)  we  will  consider  detailed  representation  of 
particular  disturbances  in  terms  of  these  eigenfunctions  and  their  spatial 
and  temporal  evolution.  We  also  plan  to  study  at  some  later  time  the 
mechanisms  (suggested  in  section  4)  by  which  these  modes  might  lead  to 
transition. 

2.  GENERAL  ANALYSIS 

2.1.  Formulation 

The  basic  flow  is  a parallel  shear  flow  (U  (y) , 0,  0),  in  a Cartesian 
coordinate  system  (x,  y,  z) . We  consider  infinitesimal  two-dimensional 
disturbances  to  this  flow.  The  stream  function  of  the  disturbance  is 
assumed  to  be  of  the  form 

^ (x,  y,  t)  = 4i  (y)  . (1) 


As  is  well  known,  with  these  assumptions  the  linearized  Navier-Stokes  equations 
reduce  to  the  Orr-Sommerfeld  equation: 


(2) 


All  variables  are  dimensionless;  the  length  scale  is  L , the  velocity  scale 
is  and  the  time  scale  . As  usual,  a is  the  wave  number,  c the 
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phase  speed,  and  R the  Reynolds  number,  L/v  , with  v the  kinematic 
viscosity.  Equation  (2)  is  to  be  solved  with  suitable  boundary  conditions 
(discussed  below) . In  the  temporal  stability  problem  a is  real  and  c is 
complex,  so  that  the  flow  is  unstable  if  Im  (c)  > 0 . For  the  spatial 
stability  problem, 

0)  = ac  (3) 

is  real  and  a is  complex;  the  flow  is  unstable  if  Im  (a)  < 0 . 

The  disturbance,  equation  (1) , is  just  a single  eigenfunction  of  the 
stability  problem.  If,  as  is  the  case  for  the  temporal  stability  problem 
of  plane  Poiseuille  flow,  the  Orr-Sommerfeld  equation  has  an  infinite  set  of 
discrete  eigenvalues  {c^^}  and  a corresponding  complete  set  of  eigenfunctions 
(y)}  > the  most  general  solution  of  the  linearized  Navier-Stokes  equations 
has,  at  fixed  R , a stream  fmction  of  the  form, 

00  00 

'f  (X,  y,  t)  = J*  ^ (a)  (x,  y,  t;  a)  da 

-*  n=l 

- f f.  \ O)  ♦„  (y;  «)  da  , (4) 

n=l 

where  the  (A^  (“) ) are  "Fourier"  amplitudes.  A physical  disturbance  in 
the  flow  will,  in  general,  excite  a mixture  of  modes,  and  this,  for  the 
case  of  temporal  stability,  corresponds  to  specifying  (x,  y,  o)  . If  any 
mode  (x,  y,  t;  a)  grows  in  time,  then  the  initial  disturbance,  or  rather 
part  of  it,  grows  in  time. 

In  order  to  generate  only  a single  mode  with  a given  wave  niomber  , 

'V  (x,  y,  o)  must  be  of  the  form 
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-IM. 


1'  (x,  y,  o)  = A 4.  (y;  o ) e^“o^  , (5) 

n o 

with  A some  arbitrary  amplitude.  It  is  difficult,  at  best,  to  conceive  of 
a physical  experiment  which  could  generate  an  initial  disturbance  exactly 
of  the  form  of  (5) . 

Actually  the  temporal  stability  problem  does  not  represent  the  controlled 
experiments  that  are  carried  out  in  shear  flows.  These  experiments  involve 
the  generation  of  a disturbance  of  fixed  frequency  at  some  fixed  x position, 
say  x^  , and  observation  of  the  spatial  growth  or  decay  of  the  disturbance, 
in  this  case 


'P  Cx,  y,  t)  = f T A (o>)  i (y;  u)  dm  . (6) 

n=l  " " 

In  principle  'll  (x^,  y,  t)  , which  is  of  the  form 

(Xq,  y.  t)  = F (y)  e^V  , (7) 

is  known,  and  so  the  {A^  (o>) } can  be  found.  If,  as  appears  to  be  the  case 
for  the  Blasius  flow,  the  higher  T-S  modes  are  highly  danqied  compared  to  the 
fundamental  T-S  mode,  they  decay  rapidly  with  increasing  x , and  some 
distance  downstream  from  the  disturbance  generator  only  the  fundamental  mode 
is  observed. 

Clearly  the  T-S  modes  are  only  a mathematical  device  for  representing 
physical  disturbances.  It  is  only  an  accident  that  the  higher  T-S  modes  are 
highly  damped  so  that  some  distance  downstream  from  a disturbance  generator 
the  pure,  fundamental  T-S  wave  can  be  observed,  even  though  a pure,  fundamental 
T-S  wave  is  not  generated  by  the  physical  disturbance  generator.  If  the 
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Orr-Sommerfeld  equation  has  eigenfunctions  other  than  the  T-S  modes, 
they  too  must  be  regarded  as  mathematical  representations  of  portions  of 
physically  realizable  disturbances. 


I 


I 


k 


For  a general  shear  flow  it  is  conceivable  that  the  Orr-Sommerfeld 
equation  might  not  have  an  infinite  set  of  discrete  eigenvalues  and  so  could 
not  possibly  have  a complete  set  of  corresponding  eigenfunctions  or  that, 
even  if  it  had  an  infinite  set  of  discrete  eigenvalues,  the  corresponding 
infinite  set  of  eigenfimctions  might  not  be  complete.  In  either  of  these 
cases  completeness  could  require  the  inclusion  of  "improper  eigenfunctions" 
(Friedman,  1956)  corresponding  to  a continuous  spectrum.  If  this  were  the 
case,  then  the  stream  function  could  not,  in  general,  be  represented  as  a 
sum  over  only  the  discrete  modes,  as  in  equation  (4) ; the  contribution  from 
the  continuum  must  be  included,  i.e. , for  the  temporal  problem 

'■  = / { E (a)  Cy;  a)  e 

-00  n 

+ B (o,  c)  (|)  (y;  a,  c)  dc  } da  , (8) 

where  the  integral  in  the  second  term  is  taken  over  the  continuum  of  c . 

This  latter  case  is  quite  common  in  other  branches  of  physics.  Consider, 

for  example,  the  nonrelativistic  Schrodinger  equation  (Landau  and  Lifschitz, 
1958)  for  an  eleccron  in  a potential.  If  the  potential  is  that  of  Hydrogen, 
it  is  easily  shown  that  there  is  an  infinite  set  of  discrete  eigenvalues  and 
corresponding  eigenfunctions  which  are  not  complete.  Completeness  requires 
the  inclusion  of  the  continuous  spectrum  and  "improper"  eigenfunctions.  If 
the  potential  is  that  of  the  Deuteron,  there  is  only  one  discrete  eigenvalue, 
?nd  completeness  requires  the  inclusion  of  the  continuum  modes  (Bethe,  1947). 
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In  both  these  cases  the  discrete  eigenvalues  are  the  energy  levels  of  the 
bound  states  of  the  electron,  and  a point  on  the  continuum  is  the  energy 
of  an  electron  which  is  scattered  by  the  potential. 


A trivial  example  may  serve  to  illustrate  the  essential  point  of  this 
argument.  Consider  the  wave  equation 


3^u  3^u 


3t2  3x2 


(9) 


If  we  look  for  solutions  of  the  form 


u (x,  t)  = f(x)  e 


ioit 


(10) 


then  f(x)  is  a solution  of 


d^f  . 9 


dx2 


+ f = 0 . 


(11) 


If  the  boundary  conditions  for  equation  (9)  are 


u (o,  t)  = u (1,  t)  = 0 , 


(12) 


then  there  is  an  infinite  set  of  discrete  eigenvalues  {uj^}  and  eigenfunctions 


(fn  (X)} 


“n  = nn  . 


(13) 


f (x)  = — sin  (nllx)  , n = 1,  2,  3,  . . . , 
n ^ 


(14) 


and  these  eigenfunctions  are  orthonormal  and  are  a complete  set.  If  we 
consider  the  infinite  domain  0 < x < » , instead  of  the  finite  domain 
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0 £ X ^ 1 , and  impose  the  boundary  conditions 


u (o,  t)  = 0 , 

(ISa) 

u (x,  t)  0 as  X -»■  » , 

(15b) 

it  is  obvious  that  there  are  a finite  number  of  discrete  eigenvalues;  in 
fact,  there  are  exactly  zero  eigenvalues.  If,  however,  the  boundary  condition 
(15b)  is  relaxed  to 


u (x,  t)  is  bounded  as  x -*•  » , 


CISC) 


then  the  spectrum  is  a continuiun, 

(i>  real  and  co  ^ 0 , 
and  the  "improper"  eigenfunctions  are 


i 

i 1 

I f(x;  (o)  * sin  tox  . (16) 

[ ^ 

I 


These  eigenfunctions  are  improper  in  the  sense  that  they  are  not  square 
integrable  in  [0,  »),  but,  as  written,  they  are  6 function  normalized,  i.e.. 


t 

I 

[ 


/ 


f(x;  (u)  f(x;  0)')  dx  = 6 (o)  - oj') 


(17) 


This  example  is,  of  course,  trivial.  The  Fourier  sine  series  on  [0,  1] 
has  been  replaced  by  the  Fourier  sine  integral  on  [0,  «).  The  fact  that  the 
continuum  eigenfunctions  are  improper  (not  square  integrable  in  [0,  “))  is  so 
familiar  as  to  be  accepted  without  comment.  The  fact  that  f(x;  oj)  is  not  a 
physically  realizable  mode  but  only  a mathematical  representation  of  a 
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portion  of  a wave  packet,  which  is  physically  realizable,  is  a familiar 
concept . 


It  will  be  shown  below  that  the  Orr-Sommerfeld  equation  for  a wide  class 
of  shear  flows  in  an  infinite  domain  always  has  a continuous  spectrum  and 
that  the  corresponding  continuum  or  "improper"  eigenfunctions  are  6 function 
normalizable.  The  first  step  in  showing  this  requires  an  examination  of  the 
boundary  conditions. 

2.2.  Boundary  Conditions  for  the  Orr-Sommerfeld  Equation 

The  Orr-Sommerfeld  equation  (2)  is  to  be  solved  with  suitable  boundary 
conditions  in  [0,  o») . First  consider  the  boundary  conditions  at  y = 0 . 

If  the  flow  is  a boundary  layer,  i.e.,  there  is  a solid  wall  at  y = 0 , 
then  the  velocity  components  (u,  v,  o)  of  the  disturbance  must  vanish  at 
y = 0 . Since 


U = V = - 

3y  ’ 3x 


(18) 


two  of  the  boundary  conditions  for  a boundary  layer  are  then 


= ^=  0 at  y = 0. 


(19) 


If  the  flow  is  an  unbounded  shear  flow  such  as  a jet,  wake  or  free  shear 
layer,  <p  can  be  written  as  a sum  of  symmetric  and  antisynmetric  modes 
(symmetry  about  y = 0 ) . If  U (y)  is  symmetric  about  y = 0 , as  is 
usual  for  a jet  or  wake,  the  symmetric  and  antisymmetric  modes  are  uncoupled 
solutions  of  equation  (2).  However,  if  U (y)  is  antisymmetric,  a free  shear 
layer  for  example,  the  symmetric  and  antisymmetric  modes  are  coupled.  For  (}i 
symmetric  about  y = 0 , the  boimdary  conditions  are 


4 


(20) 


^ = 0 


dy 


dy> 


at  y = 0 , 


while  the  boundary  conditions  for  an  antisymmetric  mode  are 


(fr  =>  =0  at  y = 0 . 

dy2 


(21) 


What  boundary  conditions  are  to  be  applied  at  infinity?  It  is  almost 
universal  to  require  that  u and  v vanish  at  infinity  and  hence  that 


’ ^ 0 as  y ->■  » 


(22) 


This  condition  insures,  of  course,  that 


/ 


(p*  (p  dy  < 00  ^ 


(23) 


with  the  star  denoting  conq)lex  conjugate,  i.e.,  (p  , is  in  L2  . This 
boundary  condition  will  not  be  used  in  this  paper.  Instead  the  weaker  condition 


<P 


dtp 

w 


are  bounded  as 


y •>  00 


(24) 


will  be  imposed. 

Clearly  if  u,  y -*■  0 as  y ■*■  <*>  they  are  also  bounded  as  y -*■  ao  . 

But,  if  there  are  modes  that  are  bounded  as  y -►  » which  do  not  satisfy 
(22)  (u,  V 0 as  y -*•  ») , then  they  are  "improper”  eigenfunctions  and  must 
be  regarded  as  a mathematical  device  for  representing  a class  of  physical 
disturbances.  It  will  be  shown  below  that  the  eigenvalues  of  the  modes  which 
satisfy  (24)  but  not  (22)  form  a continuum,  and  the  sum  over  modes  in  (4)  and 
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(6)  must  be  replaced  by  a sum  over  the  discrete  modes  and  an  integral  over  the 
continuum  as  in  equation  (8) . 

2.3,  The  Continuous  Spectrum 

2.3.1.  The  Temporal  Stability  Problem.  Equation  (2)  is  a fourth  order 
linear  differential  equation  and  therefore  has  four  linearly  independent 
solutions:  (y)  , j = 1,  2,  3,  4 . The  asynq)totic  (i.e.,  y -►  «)  form 

of  these  solutions  can  be  found  by  considering  equation  (2)  as  y -►  « . 

For  the  class  of  shear  flows  we  are  considering 

U (y)  -*■  Ui  (a  constant)  , 


U'  Cy)  - 0 , 


(25) 


and 


U"  (y)  ^ 0 , 


as  y . If  U (y)  does  not  approach  a constant  as  y , then  there 
is  no  Galilean  transformation  under  which  the  flow  has  finite  energy.  The 
constant  Uj  is  unity  for  a boundary  layer,  a wake  or  a shear  layer  and  is 
zero  for  a jet. 

Then,  as  y , equation  (2)  becomes 


(26) 


a fourth  order  differential  equation  with  constant  coefficients.  The  four 
independent  solutions  of  the  Orr-Sommerfeld  equation  are  then  asymptotic  to 
the  solutions  of  equation  (26),  i.e.. 
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'I'j  (y)  ~ , 


Xj  = - 


(28a) 


X2  = + (1^2. 


(28b) 


X3  = - a 


(28c) 


X4  = + o 


(28d) 


Q = io  R (Ui  - c)  + a' 


and  <(>2  are  the  "viscous"  solutions  while  <|»3  and  ((>4  are  "inviscid" 
solutions.  For  the  temporal  stability  problem  and  a general  complex  c , 

A 2 has  a negative  real  part  and  A2  has  a positive  real  part,  so  only  (pi 
and  <^3  satisfy  the  boundary  conditions,  equation  (24), 


<p  , bounded  as  y -»■  » . 


A linear  combination  of  and  (^>3  , the  decaying  viscous  and  inviscid 
solutions,  also  satisfy  the  more  stringent  boimdary  conditions,  equation  (22). 

The  Orr-Sommerfeld  equation,  (2),  can  then  be  solved  for  (|)2  and  (|)3  . 
Applying  the  boundary  conditions  at  y = 0 , the  eigenvalue  problem  is  reduced 
to  finding  a value  of  c for  which  a linear  combination  of  and  (j»3  will 

satisfy  two  boundary  conditioi.s,  equation  (19),  (20),  or  (21)  at  y = 0 . 

Note  that  all  these  solutions  satisfy  equation  (22)  as  well  as  equation 
(24),  i.e.,  they  are  in  L2  . However,  there  is  another  class  of  solutions 
which  satisfy  equation  (24)  but  not  equation  (22) . To  see  the  asymptotic 
(y  ->•  00)  form  of  t>ese  solutions,  assume  that  Q is  real  and  negative,  i.e.. 


Re  (Q)  = aRc^  + < 0 , 


Im  (Q)  = oR  (U^  - c^)  = 0 


From  equation  (30), 


< - o/R  , 


which  is  written,  for  convenience. 


= - (1  + k^)  a/R  , k real  and  nonzero 


Since  oR  0 , equation  (31)  gives 


= Ui 


The  trivial  solution,  a = 0 , implies  (}i  = 0 
Now  with 


c = Uj  - i (1  + k2)  a/R  , 


Q = - k2  a2  , 


Xi  = - ika  , 


(37a) 


Xo  = + ika  . 


(37b) 


Therefore,  both  viscous  solutions  (fij  and  <p2  , as  well  as  ipj  , the 
decaying  inviscid  solution,  satisfy  the  boundary  condition,  equation  (24) , as 
y -►  ® , and 
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= A <^i  + B (^2  C <t>3 


Ae 


•ikay 


* Ce-“y 


(38) 


The  independent  solutions,  (fi  , <t>2  , and  ((13  , to  equation  (2)  are 
found,  and  the  boundary  conditions  at  y = 0 , equation  (19) , (20) , or  (21) 
are  applied.  Because  there  are  three  arbitrary  constants  and  only  two  conditions 
to  be  satisfied,  the  boundary  condition  at  y = 0 can  always  be  satisfied. 

Tae  temporal  stability  problem  for  the  Orr-Sommerfeld  equation  has, 
therefore,  a continuous  spectrum  along  the  line  in  the  complex  c plane: 

c = Ui  - i (1  + k2)  a/R  , (39) 

for  arbitrary  real  positive  k . Now,  Ui  « 1 for  a (suitably  normalized) 
boundary  layer,  wake  and  shear  layer  and  1/3=0  for  a jet.  All  of  these 
continuum  modes  are  damped  (c^^  < 0);  however,  the  damping  rate, 

oc^  = - (1  + k2)  a^/R  , (40) 

is  quite  small,  at  least  for  small  k , because  o^/R  « 1 in  most  situations. 
The  continuum  eigenfunctj ins  are  the  form  of  equation  (38) , where  two  of  the 
three  constants  are  determined  from  the  boundary  conditions  at  y = 0 and  the 
third  is  arbitrary. 

When  k = 0 , 


Xj  = X2  ® 0 


(41) 


and 


' 1 » 


(42a) 
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<1,2  - y 


$ 


(42b) 


and  is  once  again  a linear  combination  of  (|ii  and  (j)3  . In  general  this 
value  of  c will  not  be  an  eigenvalue. 

2.3.2.  The  Spatial  Stability  Problan.  The  case  of  spatial  stability 
can  be  treated  in  an  exactly  similar  way.  In  this  case  m is  real  and  a 
is  complex  and  we  require  that 

Re  (Q)  = 0^2  - 0.2  - RUi  o.  < 0 , (43) 

and 

Im  CQ)  = 2o^  + RUi  o^  - Ru)  = 0 . (44) 

From  equation  (44) , 

“i ' i " - “i)  • («) 

Let  the  left-hand  side  of  the  inequality,  equation  (43)  be  set  equal  to 
- k2  r2/4  , k real  and  nonzero. 

Then,  substituting  for  o^  from  equation  (45),  it  is  found  that 

r2  - SI  } (46, 

with 

S = r2  (Ui2  + k2)/4  ; (47) 

therefore,  the  continuous  spectrum  for  the  spatial  stability  problem  lies  on 
the  curve  given  parametrically  by  equations  (45)  through  (47).  It  can  be  shown 
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that  the  positive  square  root  must  be  taken  in  equation  (46).  If  the  negative 


I 


sign  is  taken,  equation  (43)  implies  that  < 0 , which  is  impossible.  It 

can  further  be  shown  that  both  a >0  and  a.  > 0 . Since  a.  > 0 . all  of 
the  continuum  modes  are  damped  for  the  spatial  stability  problem. 

It  is  interesting  to  examine  a few  special  cases.  In  most  cases  of 
interest  oi/R  <<  1.  For  a boundary  layer,  wake  or  free  shear  layer,  Ui  = 1 , 
and  it  is  easily  shown  from  equations  (45)  through  (47)  that,  with  m/R  <<  1 
and  k -►  0 , 


to  , (48) 

-»  R [j  k2  + (to/R)2]  , (49) 


and 

c = 1 - i[j  k2  + (a)/R)2]/((o/R)  ; (50) 

and  for  k -»•  « , I 

I 

I 
I 

(51)  ! 

(52) 

(53) 

In  all  cases  c^  , the  phase  speed,  is  less  than  or  equal  to  the  free  stream 


-►  (o/k  , 

-►  kR/2  , 

c = 4(to/R)2  k'3  - i2  (to/R)  k’ 1 


Cy  ■ 0 (k  -»•  »)  , The  former  are  progressive  waves  moving  at  nearly  the  free 
stream  speed,  and  the  latter  are  standing  waves . 

Ui  = 0 for  a jet  and  it  is  easily  seen  that,  with  u/R  « 1 , 

as  k ^ 0 , 


o -*■  /o)R/2  (1  + i)  , 
c : 4I72R  (1  - i)  ; 
while  as  k ->  <»  , 


(54) 

(55) 


a ^ u)k"^  + iRk/2  , 


c = 4 (u)/R)2  k-3  - i2(u/R)k-l 


(56) 

(57) 


In  the  limit  k ^ , the  same  results  are  obtained  for  a and  c whether 

Ui  = 0 or  Ui  = 1.  The  reason  for  this  is  that  the  limiting  process  is  the  same, 
i.e. , lim  Ui/k  -►  0 . 


3.  APPLICATION  TO  SPECIFIC  SHEAR  FLOWS 
3.1  The  Blasius  Boundary  Layer 

The  flow  is  that  past  a semi-infinite  flat  plate,  and  the  velocity  profile 
U (y)  is  the  Blasius  boundary  layer  profile  (Batchelor,  1967).  The  velocity  and 
length  scales  chosen  in  this  case  are  Ug  , the  free  stream  speed,  and 

L (x)  = /vx/Ug  , (58) 

where  x is  the  distance  from  the  leading  edge.  With  this  choice  of  length  scale 
the  momentum,  6*  , and  boundary  layer,  6 , thickness  are  1.72  L (x)  and  5.04 
L (X),  respectively.  Since  U (y)  ^ 1 as  y ^ , ui  = 1 and  for  the  temporal 

stability  problem  the  continuous  spectrum  lies  along 

Cq.  - 1 - i (1  + k^)  a/R  , (59) 
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in  the  complex  wave  number  plane,  while  for  the  spatial  stability  problem  the 
continuous  spectrum  lies  along  the  curve  in  the  complex  wave  number  plane 
given  parametrically  by 


c^  = (oa*/|ap  , (60) 

Oy  = I j - 8]  j-  Vz  , (61) 

= j R (oj/a^  - 1)  , (62) 

ai.i 

8 = r2  (1  + k2)/4  . (63) 


Jordinson  (1971)  calculated  the  discrete  temporal  and  spatial  eigenvalues  of 
the  Orr-Soiranerfeld  equation  for  Blasius  flow  for  the  values:^ 

temporal  stability 

R = 580.0  , o = 0.179  , 

and 


spatial  stability 

R = 580.0  , u)  = 0.0652  , 


while  Mack  (1976)  calculated  the  discrete  eigenvalues  of  the  temporal  stability  prob- 
lem for  a number  of  values  of  a and  R including  a = 0.179  , and  R = 580.0  . 

^ The  values  given  here  differ  from  those  in  Jordinson 's  paper  because  he  used 
6*  as  the  length  scale,  instead  of  L , as  given  in  equation  (58).  Mack 
used  L as  the  length  scale. 
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Ihese  cases  lie  within  the  neutral  stability  curves  in  the  (a,  R)  plane; 
i.e. , there  is  an  unstable  T-S  mode  for  both  the  temporal  and  spatial  stability 
problem.  They  lie  near  the  top  of  the  unstable  region  about  halfway 
between  branch  I and  branch  II  (Obremski  et  al.,  1969). 

Figure  1 shows  the  coiq>lex  c plane.  The  continuous  spectra  for  both 
the  teiq)oral  problem  (a  » 0.179,  R > 580.0)  and  the  spatial  problem 
(b)  « 0.0652,  R « 580.0)  are  shown  as  well  as  the  T-S  eigenvalues  confuted  by 
Jordinson  and  Mack.^  There  is  substantial  disagreement  between  Jordinson  and 
Mack  as  to  the  number  and  location  of  the  discrete  T-S  eigenvalues.  None  of 
the  T-S  modes  found  by  Mack  lie  on  the  continuous  spectrum. 

Most  of  the  decaying  spatial  modes  found  by  Jordinson  appear  to  lie  on  a 
curve  roughly  parallel  to  the  spatial  continuous  spectrum.  It  is  difficult 
to  see  how  Jordinson  could  have  found  any  points  on  the  spatial  continuum 
since  his  boundary  condition  as  y • (equation  (3)  of  his  paper)  is 
equivalent  to  requiring  that  and  its  derivatives  vanish  as  y , and 
the  continuum  modes  do  not  approach  zero  but  only  remain  bounded  as  y . 

We  have  calculated  the  continuum  eigenfunctions  of  the  temporal  problem, 
for  the  values  of  (a,  R)  given  above,  at  a large  number  of  points  on  the 
continuous  spectrum.  The  numerical  method  is  very  similar  to  that  used  by 
Mack.  For  example,  consider  the  temporal  problem  at  a given  a and  R . 

A value  of  k is  chosen,  and  from  equation  (59)  a point  on  the  spectrum  is 
determined.  The  upper  boundary  condition  is  applied  at  some  large  y , say 
yo  , where  equation  (25)  is  satisfied  to  at  least  one  part  in  10^  . A value 

Jordinson' s computed  values  were  estimated  from  figure  1 of  his  paper,  and 
Mack’s  were  taken  from  table  la  of  his  paper. 
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of  Yo  = 9 is  sufficient  at  a = 0.179  , R = 580  ; then  at  /q  : 

♦ l * , C64) 

♦ 2 - , (65) 

and 

<(>3  = , (66) 

From  these  equations  ' (yo)  , (p."  (yg)  , and  (yg)  are  determined. 

The  Orr-Sommerfeld  equation,  (2),  is  then  numerically  integrated,  using  a 
fourth  order  Runge-Kutta  method,  to  y = 0 . A linear  combination  of  these 
three  independent  solutions  is  then  chosen  to  satisfy  the  boundary  ccmditions, 
in  this  case  equation  (19) . This  leaves  one  of  the  three  constants  undeter- 
mined. Choosing  this  constant  effectively  determines  the  normalization.  The 
normalization  we  have  chosen  is  to  set  A,  the  coefficient  of  equal  to  one. 

The  nature  of  these  eigenfimctions  can  be  seen  by  considering  them  in  the 
far  field,  i.e.,  as  y -*•  ® . For  any  k , the  continuum  stream  function  is, 
as  y , 


(D  - ^i<^Cx-ky-ct)  ^ g^ia(x+ky-ct)  ^ ^,^-cxy  ^ia(x.ct) 

= g-(l+k^)a^t/R  |gia(x-ky-t)  ^ ggia(x+ky-t)  ^ ^^-ay  ^ia(x-t)  | 

(67) 

The  amplitude  of  the  stream  function  is  decaying  exponentially  in  time  with 
a decay  time  of  R/[a^  (1  + k^)]  . Apart  from  this  decaying  amplitude,  the 
first  term  in  the  bracket  is  a progressive  wave,  with  an  amplitude  of  1, 
traveling  toward  the  solid  boundary.  The  magnitude  of  the  wave  number  is 
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o /l  + k' 


and  the  propagation  vector  lies  at  an  angle,  - 0 


with 


respect  to  the  x-axis,  where 


0 = tan“l  (k)  . (68) 

The  second  term  is  an  outgoing  wave  with  complex  amplitude  B.  The  wave 
number  is  the  same  magnitude  as  the  incoming  wave,  and  the  propagation  vector 
is  at  the  angle  + 0 with  respect  to  the  x-axis.  The  third  term  is  a 
"wall  wave"  with  complex  amplitude  C,  which  decays  outwards  from  the  wall. 

This  wave  is  propagating  parallel  to  the  wall  with  wave  number  a . 

In  the  language  of  scattering,  these  continuum  modes  consist  of  an 
incoming  wave  (amplitude  = 1 and  angle  of  incidence  = 0)  which  interacts 
with  the  boundary  layer  and  reflects  as  an  outgoing  wave  (complex  amplitude  = B 
and  angle  of  reflection  = 0) . A "wall  wave"  (compie:;  .^plitude  = C)  is 
generated  which,  in  combination  with  the  incident  and  reflected  waves,  satisfies 
the  boundary  conditions  at  y * 0 . The  complex  amplitudes  B and  C are,  of 
course,  determined  by  the  detailed  shape  of  the  velocity  profile  through  the 
Orr-Sommerfeld  equations  and  the  boundary  conditions  at  y = 0 . For  a 
fixed  a and  R , B and  C also  vary  with  k . 

The  magnitude  and  phase  of  the  complex  amplitudes  B and  C are  shown  in 
figure  2 as  a function  of  0 , the  angle  of  incidence  of  the  incoming  wave, 
for  the  temporal  stability  continuum  at  a = 0.179  , R = 580.0  . The 
magnitudes  of  B and  C are  plotted  in  figure  2a,  and  the  phases  of  B and  C are 
plotted  in  figure  2b.  Both  the  magnitudes  and  phases  show  somewhat  complicated 
variation  with  0 . 

As  0 -*•  0 , i.e.,  at  grazing  incidence  [b]  ->■  1 , phase  (B)  -►  n , 

lc|  -*•  0 , and  phase  (C)  - 99°  . As  0 increases,  |b|  and  phase  (B) 
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both  decrease  slowly  up  to  0 =65®  . In  this  same  range  |c|  increases  and 
the  phase  (C)  decreases.  The  maximum  of  |c|  occurs  at  0 = 65®  , where 
|c|  = 3.75  and  phase  (C)  ~ -311/2  . For  larger  0 , [cj  decreases  very 
rapidly,  and  both  phase  (B)  and  phase  (C)  undergo  rapid  variation.  The 
variations  of  phase  (B)  and  phase  (C)  are  only  shown  up  to  0 - 85®  . The 
variation  of  the  phase  with  0 is  so  rapid  for  larger  angles  of  incidence  that 
it  is  not  possible  to  show  it  on  this  figure. 

The  real  and  imaginary  parts  of  are  shown  in  figure  3 as  a function 
of  y for  the  case  a = 0.179  , R = 580.0  , and  k = 2.0  (0  = 63.43®). 

It  can  be  seen  that  the  disturbance  is  very  small  within  the  lower  part  of 
the  boundary  layer  (the  top  of  the  boimdary  layer  is  at  y = 5 . 04) . This  was 
found  to  be  true  for  all  0 . These  continuum  modes  are  essentially  free 
stream  modes  and  do  not  penetrate  very  far  into  the  boundary  layer. 

Contours  of  the  stream  function,  , in  the  far  field  for  this  same 
case  are  shown  in  figure  4.  The  region  shown  is  0 ^ ax  £ 4II  and 
5n  < ay  ^ 7n  . It  can  be  seen  that  these  continuum  modes  are  a doiobly 
periodic  airay  of  vortices  in  the  free  stream.  The  ratio  of  the  wavelength 
in  the  x-direction  to  that  in  the  y-direction  is  k , in  this  case  2.0. 

Because  |a|  ^ |B|  , the  axes  are  tilted  and  the  vortices  are  distorted. 

The  vortex  array  is  moving  with  the  free  stream,  i.e.,  has  a phase  speed  of 
1.0,  and  the  vortex  strength  is  decaying  in  time  as  exp  [-  (1  + k^)  a^t/R]  . 

It  is  readily  seen  that  all  the  continuum  modes  are  of  the  same  general  form 
in  the  far  stream.  The  aspect  ratio  of  the  vortices  and  the  distortion 
depend  on  k and  B (k)  . 

Some  related  results  have  been  obtained  by  Rogler  (1975)  and  Rogler 
and  Reshotko  (1975)  who  considered  the  response  of  a boundary  layer  to  an 
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array  of  nondecaying  free  stream  vortices.  This  problem  is  not,  as 
emphasized  by  Rogler  and  Reshotko,  an  eigenvalue  problem,  because  the 
form  of  the  vortex  array  is  imposed  in  the  far  field,  as  well  as  the 
requirement  that  the  vorticity  be  constant  in  time.  Applying  a perturbation 
analysis  to  the  Navier-Stokes  equations,  Rogler  and  Reshotko  obtain,  formally, 
an  inhomogeneous  Orr-Sommerfeld  equation  with  a forcing  function  depending  on 
the  mean  velocity  field  and  the  vortex  array.  It  should  be  en^hasized  that 
the  vortex  arrays  considered  by  Rogler  and  Reshotko  are  not  the  continuum 
eigenfunctions  of  the  Orr-Sommerfeld  equation.  The  connection  between  these 
continuum  modes  and  the  free  stream  vortices  used  by  Rogler  and  Reshotko 
will  be  discussed  in  a later  paper  (Part  II). 

3.2.  The  Two-dimensional  (Plane)  Laminar  Jet 

As  a second  example,  the  continuum  modi  s of  the  two-dimensional  laminar 
jet  will  be  discussed.  The  velocity  of  this  jet  is,  in  dimensionless  units. 


U (y)  = 1 - tanh^  (y)  . 


(69) 


The  length  sccle,  L (x)  , and  velocity  scale,  Uq  (x)  , are 


L (x)  = [48  x2/K]^  , 


(70) 


and 


Ug  (X)  = [3K2/32  V x]l^  , 


(71) 


where,  the  cap  (""”)  denoting  a dimensional  quantity. 


4*^  A A 

K = / U (y)  dy 


(72) 


! I 
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is  the  momentum  flux  per  unit  density,  per  unit  width  of  the  jet  (Schlichting, 
1951). 

It  is  clear  that  the  parallel  shear  flow  approximation  is  poorer  for  a 
jet  (L  (x)  - x^^)  than  for  a boundary  layer  (L  (x)  - x^) , but  we  will  not 
consider  here  the  modification  to  the  continuum  necessary  in  the  nonparallel 
flow  approximation  (but  see  below,  section  4) . 

The  free  stream  speed  is  zero  for  the  jet,  and  so  the  continuous  spectrum 
for  the  temporal  stability  problem  lies  along  the  line 


c„  = - i (1  + k^)  a/R  . 


(73) 


As  discussed  above,  this  symmetric  jet  has  two  continuum  modes,  the 
symmetric  and  antisymmetric  modes,  for  each  point  on  the  continuum.  In  both 
cases  the  modes  are  decaying,  standing  waves. 

The  continuum  of  the  spatial  stability  problem  for  the  jet  lies  along 
the  curve 


= (k  R/2)  I ^ [/iTi-  16  (u)/k^R)^  - 1]  } ^ . 
“i  " {t  [^1  + 16  (w/k^R)2  - 1]  . 


(74) 


(75) 


In  contrast  to  the  temporal  stability  problem,  the  continuum  modes  of  the 
jet  for  the  spatial  stability  problem  are  spatially  decaying,  traveling  waves 
with  a wave  speed. 


(o/otj.  = (2o)/k  R)  [/I  + 16  (u)/k^R)2  - i]|  -1/2  . 


(76) 
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We  have  calculated  the  symmetric  and  antisymmetric  continuum  modes  of 
the  laminar  jet  for  the  temporal  stability  problem  at  a = 1.0  , R = 50.0  , 
and  0 < k ^ 100.0  , (0®  < 0 s 89.4®)  . The  numerical  method  is  identical 

to  that  used  with  the  boundary  layer. 

Although  the  form  of  the  eigenfunction,  (ji  (y)  , is  different  for  the 
symmetric  and  antisymmetric  modes,  the  complex  amplitudes  B and  C at  a given 
k are  nearly  identical.  The  amplitude  and  phase  of  B and  C for  both  the 
symmetric  and  antisymmetric  modes  are  plotted  versus  6 in  figure  5.  The 
differences  between  the  results  for  the  symmetric  and  antisymmetric  are  so 
small  that  they  cannot  be  shown  on  this  figure. 

The  variations  of  the  amplitudes  and  phases  of  B and  C are  much  smoother 
for  the  jet  than  for  the  boundary  layer.  As  0 -►  0 , amp  (B)  -►  1.0  and 
amp  (C)  -►  0 , while  phase  (B)  -►  180®  and  phase  (C)  -►  - 45®  . The 
amplitiides  of  both  B and  C rise  monotonically  with  increasing  0 to  a peak 
response  at  0 = 79.5®  (k  = 5.15)  and  then  decreases  rapidly  thereafter.  The 
phases  of  B and  C decrease  monotonically  with  increasing  0 , slowly  for 
0 < 79.5®  and  then  very  rapidly.  The  phases  are  not  shown  for  0 > 85®  because 
the  extremely  rapid  variation  with  0 in  this  region  makes  it  impractical. 

The  maximum  amplitudes  of  B and  C are  very  large  for  the  jet  as  compared 
to  the  boundary  layer  when  it  is  noted  that  in  both  cases  the  normalization 
is  the  same,  i.e.,  A = 1 for  all  0 . In  the  case  of  the  jet,  the  maximum 
amplitude  of  B is  about  4850  (1.15  for  a botmdary  layer),  and  the  maximum 
amplitude  of  C is  about  2300  (3.75  for  a boundary  layer).  This  was  to  be 
expected  in  view  of  the  very  low  critical  Reynolds  number  for  a jet  (of 
order  10)  as  compared  to  the  critical  Reynolds  number  for  the  Blasius  boundary 
layer  (about  300  with  this  length  scale) . 
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This  same  sensitivity  of  the  jet  to  the  free  stream  disturbances  can  be  \ 

seen  in  figures  6 and  7 wh^re  the  Re  (^)  and  Im  are  plotted  versus  y ] 

for  the  symmetric  and  antisymmetric  modes  with  a^l.O,  R>  50.0,  and  j 

k s 2.0  (6  > 63.43°)*  It  can  be  seen  that  in  both  cases  the  disturbance 

f. 

penetrates  much  further  into  the  jet  than  was  the  case  for  the  boundary  layer. 

If  the  "edge"  of  the  undisturbed  jet  is  defined  as  the  point  where  U Cy)  ~ 0.001  , 
then  the  jet  occiq>ies  the  region  |y|  S3,  and  it  is  seen  from  figures  6 and 
7 that  there  is  quite  a substantial  disturbance  in  this  region.  If  the  an;)li- 
tude  of  the  disturbance  were  large  enough  the  apparent  "edge"  of  the  jet  could 
be  moved  by  the  disturbance. 

Figure  8 shows  the  far  field  stream  function  for  the  symmetric  mode  of 
the  jet  at  a « 1.0  , K » 50.0  , and  k > 2.0  . The  region  shown  is 
0 <,  ouc  £ 8n  and  lOII  ^ ay  ^ 14n  . As  in  the  case  of  the  Blasius  boundary 
layer,  the  far  field  is  a doubly  periodic  array  of  vortices,  but  in  contrast 
to  the  boundary  layer  the  vortices  are  highly  distorted  due  to  the  fact  that 
|b|  is  considerably  different  from  |a|  > 1 . At  this  distance,  y > 30  , 
from  the  jet  axis  the  trapped  mode  is  negligible.  Ihe  far  field  stream 
function  for  the  antisymmetric  mode  is  essentially  identical  to  that  of  the 
symmetric  mode  and  is  not  shown. 

4.  THE  CONTINUUM  MODES  AND  TRANSITION 

All  of  the  continuum  modes,  both  teiq>oral  and  spatial,  for  the  linearized 
stability  problem  of  an  unbounded  parallel  shear  flow  are  damped;  the  flow  is 
linearly  stable  with  respect  to  these  modes.  It  is,  of  course,  possible 
that  if  the  free  stream  vorticity  is  sufficiently  intense  the  nonlinear 

j 

problem  would  yield  growing  solutions,  i.e.,  the  flow  would  be  unstable  to 
these  modes  and  transition  would  result. 

j 
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We  believe,  however,  that  there  are  three  possible  mechanisms  by  which 


small  amplitude  continuum  modes  might  lead  to  transition.  The  first 
mechanism  involves  a small  quasi-steady  distortion  of  the  mean  profile  by 
the  continuum  modes;  the  second  mechanism  requires  the  consideration  of  the 
fully  time -dependent  flow,  consisting  of  the  mean  flow  plus  the  slowly 
decaying  continuum  modes;  while  the  third  mechanism  involves  consideration 
of  the  nonparallel  flow  corrections  to  the  continuum  spectrum. 

The  usual  picture  of  transition  via  T-S  waves  in  a boundary  layer  consists 
of:  the  appearance  of  a small  amplitude,  unstable  T-S  wave  of  unspecified  origin; 

its  growth  through  the  region  of  validity  of  the  linear  theory;  the  appearance 
of  small  nonlinear  effects,  i.e.,  higher  harmonics  and,  more  importantly,  a 
small  three-dimensional  distortion  of  the  mean  profile  near  the  top  of 
the  boundary  layer;  and  finally  the  appearance  of  short  wavelength,  rapidly 
growing  T-S  waves  which  lead  directly  to  the  appearance  of  a turbulent 
burst.  It  is  believed  that  the  penultimate  stage,  the  distortion  of  the  mean 
profile  near  the  top  of  the  boundary  layer,  is  crucial.  This  distortion, 
in  the  region  where  U”  is  small,  causes  U"  to  change  sign,  and  the 
inflection  in  the  mean  profile  permits  the  appearance  of  small  wavelength, 
highly  unstable  T-S  waves.  These,  presumably,  become  nonlinear  rather 
quickly  and  cause  a further  distortion  of  the  mean  profile,  and  so  on  to 
transition.  This  is,  in  effect,  a modified  Landau  model. 

Now  consider  the  effect  of  a small  patch  of  vorticity  in  the  free  stream. 
Assuming,  as  we  will  show  in  Part  II,  that  we  can  represent  this  by  an 
integral  over  the  continuum,  all  wavelength  continuum  modes  are  present  and 
exciting  the  boundary  layer.  Among  these  are  small  a , small  k continuum 
modes,  i.e.,  lightly  danced,  long  wavelength  and,  hence,  low  frequency  modes. 
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In  a quasi-steady  model  these  produce  a small,  quasi-steady  distortion  of 
the  mean  profile  in  the  outer  regions  where  U"  is  small;  therefore,  we 
might  expect  that  the  boundary  layer  has  reached  the  critical  stage  when 
short  wavelength,  highly  unstable  T-S  waves  can  appear  without  ever  experiencing 
a long  wavelength,  slowly  growing  T-S  wave.  In  effect,  this  mechanism  would 
bypass  the  linear  T-S  stage  and  appear  to  lead  directly  to  the  late  stages 
of  transition. 

The  second  mechanism  involves  consideration  of  the  fully  time-dependent 
flow,  consisting  of  a superposition  of  the  mean  flow  and  the  small  amplitude 
continuum  modes.  Again  consider  the  effect  of  a small  patch  of  free  stream 
vorticity.  Among  the  modes  necessary  to  represent  this  patch  there  will 
be  small  a modes,  i.e.,  long  wavelength,  lightly  damped  waves.  Let  us 
consider  the  stability  of  this  time-dependent  flow:  the  mean  flow  plus 

the  lightly  damped  modes.  Relatively  little  is  known  about  the  stability  of 
time-dependent  shear  flows  (see  Davis,  1976  for  a recent  review),  but  it  is 
known  that  under  certain  circumstances  these  flows  can  have  resonances. 

That  is,  if  the  "matching"  is  proper,  the  time-dependent  flow  can  be  unstable 
to  T-S  like  modes  which  extract  energy  from  the  mean  flow  at  a rapid  rate;  they 
are  highly  unstable.  The  unstable  mode  tends  to  lie  near  a higher  mode  of 
the  steady  Orr-Sommerfeld  equation  (Grosch  and  Salwen,  1969).  Again,  if  this 
were  to  occur,  the  linear  T-S  region  of  instability  would  be  bypassed  and 
it  would  appear,  in  an  experiment,  that  the  short  wavelength,  highly  unstable 
waves  had  been  directly  excited. 

Finally,  let  us  consider  the  possible  effects  of  nonparallel  flow 
corrections  on  these  continuum  modes.  Saric  and  Nayfeh  (1975)  have  developed 
a general  theory  which  permits  the  calculation  of  the  correction  to  the 
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spectrum  of  the  Orr-Sommerfeld  equation  due  to  the  nonparallel  nature  of 
any  nearly  parallel  shear  flow.  They  have  applied  this  theory  to  the 
stability  of  the  Falkner-Skan  flows.  In  all  cases  which  they  considered, 
the  nonparallel  flow  corrections  resulted  in  a decrease  in  the  stability  of 
the  boimdary  layer  as  compared  to  the  parallel  flow  calculation  at  fixed 
(real)  wave  number  a and  Reynolds  number  R . That  is,  at  a fixed  o , 
the  flow  was  predicted  to  be  unstable  at  a lower  R than  that  predicted 
by  the  parallel  shear  flow  approximation;  and,  at  a fixed  R , the  band  of 
unstable  wave  numbers  was  wider  than  that  predicted  from  the  usual  parallel 
flow  approximation. 

The  continuum  modes,  particularly  at  the  long  wavelength  end  of  the 
spectrum,  are  very  lightly  damped.  It  does  not  appear  to  be  unreasonable 
that  nonparallel  flow  corrections  may  make  a portion  of  the  continuous 
spectrum  unstable. 

All  of  these  mechanisms  are,  at  this  time,  only  hypotheses.  We  know 
of  no  direct  evidence  for  or  against  any  of  them,  other  than  the  observations 
that  free  stream  turbulence  can  trigger  transition  without  the  appearance  of 
the  usual  T-S  waves.  We  believe,  however,  that  all  of  these  mechanisms  are 
sufficiently-  plausible  that  experimental  and  theoretical  tests  of  them  is 
warranted.  We  hope  to  carry  out  a theoretical  investigation  of  these 
proposed  transition  mechanisms;  we  are  currently  investigating  the  nonparallel 
flow  corrections. 

5.  SUMMARY 

We  have  shown  that  the  Orr-Sommerfeld  equation  (for  both  the  temporal  and 
spatial  problems)  has  a continuous  spectrum  for  any  mean  flow  which  is  an 
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unbounded  shear  flow  and  has  finite  energy  under  some  Galilean  trans- 
formation. 


Formulae  for  the  location  of  the  continuum  in  the  complex  wave  speed 
plane  have  been  given.  These  results  have  been  applied  to  two  specific  flows: 
the  Blasius  boundary  layer  and  the  laminar  jet.  For  both  of  these  flows 
the  continuous  spectra  have  been  given.  The  continuum  eigenfunctions  of 
the  temporal  problem  have  been  calculated  for  both  flows  for  a single  wave 
number  and  Reynolds  number.  It  has  been  shown  that  the  eigenfunctions  in  the 
free  stream  are  a doubly  periodic  array  of  vortices  and  that  these 
eigenfunctions  can  be  thought  of  as  an  incident  wave,  a reflected  wave  and 
a "wall"  wave.  Holding  the  amplitude  of  the  xncident  wave  fixed,  the 
amplitudes  and  phases  of  the  reflected  and  "wall"  waves  vary  with  the  angle 
of  incidence,  and  there  is  an  optimum  angle  of  incidence  which  maximizes 
the  amplitude  of  the  wall  wave.  It  was  shown  that  the  continuum  modes  do 
not  penetrate  deeply  into  the  region  of  large  shear  in  the  main  flow  of  a 
boundary  layer  although  they  penetrate  further  and  have  much  larger 
amplitudes  in  the  jet  than  in  the  boundary  layer.  Finally,  three  mechanisms 
were  proposed  by  which  these  continuum  modes  could  cause  transition  in  a 
shear  flow  while  bypassing  the  usual  linear  T-S  stage.  These  mechanisms  are, 
at  present,  hypotheses;  we  hope  to  investigate  these  hypotheses  theoretically 
as  well  as  to  investigate  the  relation  of  these  eigenfunctions  to  the 
forced  free  stream  disturbances  studied  by  others.  We  hope  that  experiment- 
alists will  also  look  for  these  modes. 
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temporal  spectrum  of  the  Orr-Sommerfeld  equation  for  Blaslus  flow 'at  R » 580.0 
(spatial)  or  o = 0.179  (temporal).  The  points  of  the  discrete  spectrum  at  this 
a , as  computed  by  Jordinson  (1971),  are  denoted  by  0-temporal  inodes,  and 
, while  the  results  of  Mack's  calculations  (temporal  modes)  are  denoted  by  a . 
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Figure  2 


. The  magnitudes  and  phases  of  the  complex  amplitudes  B 
and  C for  the  continuum  modes  of  the  Blasius  boundary 
layer  at  a = 0.179,  R = 580. 


Figure  3.  The  real  (y)  (solid  line)  and  imaginary  (y)  (dashed 

line) , parts  of  (y)  versus  y of  the  continuum  mode  of  the 
Blasius  boundary  layer  at  a * 0.179,  R = 580,  and  k = 2. 

(0  » 63.43*) . 
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Figure  5.  The  magnitude  and  phase  of  the  con^lex  amplitudes  B and 
C for  the  continuum  modes  of  the  two-dimensional  laminar 
jet  at  o » 1 and  R = 50. 


Figure  6.  The  real  (y)  (solid  line)  and  the  imaginary  (p.  (y)  (dashed 
line),  parts  of  <t>  for  the  symmetric  temporal  stability  mode  of 
a two-dimensional  laminar  jet.  Here  o = 1,  R = 50,  and  k = 2. 

(0  = 63.43°). 


Figure  8.  The  far  field  stream  function  of  the  symmetric  continuum  mode  of  the  laminar  jet  at 
and  k = 2.  (6  = 63.43"). 
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